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ABSTRACT 

We consider the possibility that an infalling bulge core or stellar cluster could form an eccentric disk, 
following tidal disruption by a massive black hole in the center of a galaxy. As a function of central black 
hole mass, we constrain the core radii and central densities of cluster progenitors capable of becoming 
nearly Keplerian disks which can support lopsided slow modes. We find that progenitor stellar clusters 
with core radii less than a pc and densities above a few times lO 5 M0pc -3 are likely eccentric disk 
progenitors near a massive black hole of mass 10 7 - lO 8 Af . Lower density and larger progenitor cores 
are capable of causing eccentric stellar disks near more massive black holes. Our constraints on the 
progenitor cores are consistent with existing N-body simulations, which in one case has produced an 
eccentric disk. For M31 and NGC 4486B, the estimated progenitor cluster cores are dense and compact 
compared to Galactic globular clusters, however the cores of nearby galaxies such as M33, M32 and M31 
itself are in the right regime. If galaxy mergers can create eccentric disks, then they would be a natural 
consequence of hierarchical galaxy formation. 



1. INTRODUCTION 

Ever since the double-peaked nucleus of M31 was clearly 
resolved by Hubble Space Telescope (HST) imaging (Lauer 
et al. 1993), the morphology of this system has been a chal- 
lenge to explain. The most successful kinematic model 
is the eccentric stellar disk proposed by Tremaine (1995) 
where stars in apsidally aligned elliptical orbits about 
a massive black hole result in an over-density at their 
apoapses, thus accounting for the brighter peak PI. The 
black hole itself resides near the fainter, more centrally 
located peak denoted P2. Though this model has been 
successful at matching the observed velocity and luminos- 
ity distribution (Kormendy & Bender 1999; Statler et al. 
1999; Bacon et al. 2001; Sambhus & Sridhar 2002), eccen- 
tric disk formation remains a mystery. 

Touma (2002) showed that a small fraction of counter 
rotating stars, possibly originating from an disrupted glob- 
ular cluster on a retrograde orbit, could cause a more 
massive pre-existing stellar disk to develop a lopsided or 
m = 1 instability. The self-consistent kinematic modeling 
of Sambhus & Sridhar (2002), which requires a small per- 
centage of counter rotating stars, supports this scenario. 
Bacon et al. (2001) proposed that a collision from a pass- 
ing molecular cloud or globular cluster could knock a pre- 
existing stellar disk off-center, resulting in a long lived, 
precessing mode. The N-body simulations of Jacobs & 
Sellwood (1999); Taga (2002) have illustrated that lop- 
sided stellar disks could be long-lived. 

The scenarios discussed by Sambhus & Sridhar (2002); 
Touma (2002); Jacobs & Sellwood (1999); Taga (2002); 
Bacon et al. (2001), begin with an initially axisymmetric 
stellar disk, which then becomes lopsided either because of 
a violent event (such as a collision with a globular cluster) 
or the growth of an instability. However Bekki (2000a) 
proposed that M31's eccentric disk could have been a re- 
sult of a single disruption event. He succeeded at produc- 
ing an eccentric stellar disk by disrupting a globular cluster 
near a massive black hole in an N-body simulation, though 



an SPH simulation of the disruption of a massive gaseous 
cloud did not yield an eccentric disk (Bekki 2000b). The 
massive black hole in in M31 is moderate with a mass 
3 - 7 x 10 7 M Q (Kormendy & Bender 1999; Bacon et al. 
2001). The eccentric disk itself is nearly as massive as the 
black hole, - 3 x 10 7 M Q (Peng 2002; Sambhus & Srid- 
har 2002; Bacon et al. 2001), presenting a problem for 
the globular cluster disruption scenario proposed by Bekki 
(2000a) which assumed that the globular cluster was typ- 
ical of Galactic globular clusters and of order a million 
solar masses. 

Here we reconsider Bekki's proposal, that a single dis- 
ruption event could have produced the eccentric disk in 
M31. The disruption event is likely to be complex so we 
focus on what final disks are capable of supporting lopsided 
a slow modes, drawing from the recent work of Tremaine 
(2001) who developed a formalism for the purpose of pre- 
dicting the precession rates and eccentricities of discrete 
modes for low mass, nearly Keplerian disks. 

We place limits on the density and core radius of progen- 
itor clusters for eccentric disks near massive black holes. 
The resulting diagram is useful toward predicting what 
types of galaxies are most likely to harbor eccentric disks, 
for predicting the probability that eccentric disk forma- 
tion events occur and for guiding the initial conditions of 
N-body simulations which may determine if and how they 
form. 

2. TIDAL DISRUPTION OF A CLUSTER 

The observed correlation between mass of a black hole 
and the bulge dispersion (Gebhardt et al. 2000; Ferrarese 
& Merritt 2001) allows us to relate the black hole to the 
bulge of its host galaxy. 

log(M bh /M Q ) = (8.13±0.06)+(4.02±0.32) log 



,200km s" 1 
(1) 

where ex* is the stellar bulge dispersion, and we adopt con- 
stants in the relation given by Tremaine et al. (2002). For 
a discussion on the systematics of measuring the velocity 
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dispersion within r e , the effective or half light radius, see 
Tremaine et al. (2002). 

The transition radius (sometimes called the sphere of in- 
fluence) where the gravity from the black hole takes over 
from the bulge is 
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Using the above relation between the bulge and the black 
hole mass, 

We now consider disruption of a stellar cluster or molec- 
ular cloud by a massive black hole. The tidal force from 
the black hole and bulge is 
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at distance r from the galaxy nucleus and distance s from 
the cluster nucleus. We set this equal to the gravitational 
force from the cluster on itself to determine how the cluster 
disrupts. 

We describe the nucleus of a cluster with two parame- 
ters, the core radius, ro, and the density within this radius, 
po . The cores of isothermal spheres and King models are 
characterized with with these two parameters (e.g., Binney 
& Tremaine 1987). Within the core radius of the cluster, 
we assume that the cluster has nearly constant density, so 
its gravitational potential is proportional to s. We set the 
tidal force from the galaxy equal to that from the clus- 
ter and determine the distance r d , (the disruption radius) 
when the two are equal. The distance from the galaxy nu- 
cleus where the entire core of the cluster disrupts is then 
set only by the cluster core density, po- 
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Inside the transition radius the disruption radius 
r d ~ 1 | = 4.6pc [ 
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Outside the transition radius 
r d ~ 77^ = 4.7pc 
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We expect that it is difficult to make a slowly process- 
ing eccentric disk in a region where the black hole does 
not dominate the gravitational field. Outside the sphere of 
influence, the potential is strongly non-Keplerian. The re- 
quirement that the cluster be disrupted (survive to) within 
the transition radius is r d < r t . Using the previous equa- 
tion, the definition for r t , and the relationship between 
the bulge stellar dispersion and the black hole mass this 
constraint gives the requirement 

^xlO^pc-^p (,) 

From this we see that the less massive the black hole, the 
more dense a cluster would be required to form an eccentric 
disk. This condition is our first constraint on progenitor 
clusters for eccentric disks. 



At radii larger than the disruption radius, the outer 
parts of the stellar cluster or galaxy core would be 
stripped. We expect that the radial distribution of the 
resulting disk will depend upon the concentration of the 
cluster or radial profile of the infalling galaxy bulge or core. 
Since the densest part of the cluster is disrupted at r d , we 
expect that this part is most likely to be involved in any 
large eccentricity amplitude variation. So we can estimate 
the mass of the active disk formed at disruption using the 
core radius of the cluster, ro, resulting in an active disk 
of mass equal to the mass within the core of the cluster, 
M c ~ pqTq. The width of this active disk would be of the 
order of the core radius of the cluster. 

Tremaine (2001); Lee & Goodman (1999) assume that 
the mass of the disk is less than that of the black hole. If 
we consider the active part of the disk, this is even true in 
M31 (the center of mass is closest the black hole). It's dif- 
ficult to imagine a stable situation when the lopsided part 
of the disk is as massive as the black hole so we arrive at 
another constraint on our progenitor cluster, 

M c < M bh . (9) 

It is unlikely that a thin disk will result following disrup- 
tion unless the disk Mach number, M(r) > 1 (Tremaine 
2001), 
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where il(r d ) is the angular rotation rate at the disruption 
radius r d , and a c is the central cluster dispersion. 

Following disruption of the cluster, the energy spread 
of the particles is set by the cluster core radius, ro, so we 
expect that the resulting width of the disk ~ r . The an- 
gular momentum spread at disruption is also set by the 
size of the cluster core, ~ — where Jn is the angular 
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momentum of the bulk motion of the cluster at disruption. 
Consequently following disruption we expect a disk radial 
velocity dispersion at least 
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This scaling is similar to that outlined by Johnston et al. 
(2001) for the disruption of dwarf galaxies. Alternately we 
can use the relation between central dispersion, core den- 
sity and radius for a King model {a\ = 47rG ( o r ! /9). Us- 
ing either estimate for a c and inserting into equation(lO), 
assum e that we are within the transition radius so that 
we estimate 
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which is equivalent to requiring M b h > M d - 

In summary we can constrain properties of an incoming 
cluster which subsequently forms an eccentric disk: 

1. r d < r t . The disruption radius should be less than 
the transition radius. 

2. ro < rt. The cluster core radius should be less than 
the transition radius. 

3. M c < M b h- The cluster core mass should be less 
than the black hole mass. This constraint is equiv- 
alent to requiring ro < r d or M. > 1. 
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We plot these three constraints as lines in Figure 1 for 
three different black hole masses, and as a function of clus- 
ter core radius, r , and density, p . The allowed region of 
parameter space for a disrupted cluster to form an eccen- 
tric disk is the upper left hand corner, corresponding to 
core radii above about a pc, and core densities above about 
lO 5 M pc- 3 . 

Also shown in Figure 1 are measured core radii and cen- 
tral densities (assuming a mass to light ratio of M/Ly = 5 
in solar units) for the Milky Way globular clusters based 
on the quantities compiled by Harris (1996). 1 

We also include on Figure 1, the massive globular clus- 
ter Gl (based on core properties estimated by Meylan et 
al. 2001) and the cores of M31, M32 and M33, based on 
density profiles (or limits in the case of M32) estimated 
by Lauer et al. (1998). Limits on the density and core ra- 
dius are placed on the figure 1 for the two densest galaxies 
(other than M31, M32, and M33) studied by Faber et al. 
(1997). Higher angular resolution observations would be 
needed to find out if more distant galaxies can have cores 
with the high stellar densities of M31, M32 and M33. 

We see from Figure 1 that Galactic globular clusters are 
capable of forming eccentric disks if they are extremely 
dense and have compact cores, and they are more likely to 
do so in galaxy nuclei containing larger black holes. How- 
ever, to provide a plausible progenitor for an eccentric disk, 
a cluster must not only have a dense and compact core but 
be massive enough to account for the mass of the entire 
disk. Galactic globular clusters are not massive enough to 
account for the disk in M31. 

2.1. Bekki's simulations 

We examine the properties of the simulated disruptions 
done by Bckki (2000a, b) to see if they are consistent with 
our proposed progenitor cluster limits. In an N-body sim- 
ulation, Bekki (2000a) succeeded in producing a stellar 
eccentric disk by disrupting a globular cluster of mass 0.1 
times the black hole mass, for Mbh — 1O 7 M . The bulge 
he adopted was described by an NFW profile with a scale 
radius of lOOpc. For this adopted bulge profile we calcu- 
late that the mass enclosed within lOpc was less than a 
million Mq so that the entire disruption took place within 
what we have called the transition radius, where the grav- 
ity from the black hole takes over from the bulge (equation 
2). His simulated cluster disrupted at a radius of r& ~ 3pc 
from the black hole, and since the cluster core was much 
smaller than this, r§ ~ 0.4pc, the stellar system was ca- 
pable of forming an eccentric disk. We conclude that the 
success of his simulation in producing the eccentric disk is 
consistent with the crude limits we have placed on progen- 
itor cluster properties. 

The SPH simulation by Bekki (2000b), however, pro- 
duced a circular disk, not an eccentric one. The gas cloud 
in this simulation was of similar mass to that of the stellar 
cluster, (Bckki 2000a) but had a significantly larger scale 
length of lOpc. The cluster disrupted at ~ lOpc which 
is approximately the same as the core radius of the gas 
cloud. We suspect that a gas cloud with a smaller core 
radius would have been capable of producing an eccentric 
disk. Again, we find that the failure of this simulation in 
producing the eccentric disk is consistent with the crude 

1 We used the updated catalog available at http://physun.physics.rii 



limits we have placed on progenitor cluster properties. 

2.2. M31 's eccentric disk and possible progenitor cluster 
or galaxy core properties 

The two peaks observed in the nucleus of M31 are sep- 
arated by 1.87 pc (Lauer et al. 1993) and the black hole 
mass is 3 x 10 7 M Q (Kormendy & Bender 1999), so we 
expect a transition radius of rt ~ 3.1 pc which is, as we 
expect, outside the location of the eccentric part of the 
disk. The width of the eccentric part of the disk is about 
0.2" or 0.7pc so we would estimate that if the disk formed 
from a disrupted cluster that its core radius would have 
been smaller. The total mass in the most eccentric part of 
the disk is about 10 6 M Q (Bacon et al. 2001) so we would 
require a progenitor cluster of r ~ 0.7pc and po ~ a few 
times 10 6 M Q pc~ 3 . For p ~4x 10 6 M Q pc -3 the disrup- 
tion radius (estimated by equation 6) is consistent with 
the location of the eccentric disk itself. 

To provide a plausible progenitor, a cluster must not 
only have a dense and compact core but be massive enough 
to account for the mass of the entire disk. The extended, 
nearly circular part of the M31 stellar disk dominates its 
mass and has been estimated to be 3 x 10 7 M Q (Peng 
2002; Sambhus & Sridhar 2002; Bacon et al. 2001). This 
is massive compared to Galactic globular clusters, pre- 
senting a problem for the globular cluster disruption sce- 
nario proposed by Bekki (2000a). However extremely mas- 
sive extra-galactic clusters have been identified. For com- 
parison, we show in Figure 1 the location of the mas- 
sive globular cluster (or perhaps core of a dwarf galaxy) 
which is 40kpc from the nucleus of M31. This cluster 
has a mass of 7 - 17 x 1O 6 M (Meylan et al. 2001), a 
core radius r ~ 0.14" = 0.52pc, and a core density 
po ~ 4.7 x lO 5 M pc -3 (Meylan et al. 2001). Gl is mas- 
sive enough that it could have accounted for the mass in 
the extended (and almost circular) part of M31's eccentric 
disk, though its tidal radius is large, ~ 200pc, which im- 
plies that much of the cluster would have disrupted well 
outside the location of the eccentric disk in M31. 

The estimated progenitor core properties required to 
form the M31 eccentric disk are surprising if we compare 
them to Milky Way globular clusters and nearly matched 
by the exotic and extremely massive extra-galactic cluster 
Gl. However, they are similar to what has been measured 
in galaxy nearby cores, such as M33, M32 and M31 itself, 
which we have also placed on Figure 1. We see from this 
figure that galaxy cores can have the high densities and 
compactness to be eccentric disk progenitor candidates. 
Moreover, they are usually more massive than globular 
clusters and so can more easily account for the mass in 
M31's eccentric disk. 

2.3. NGC 44$6' s eccentric disk and its progenitor 

The other well known example with a double nucleus is 
NGC 4486B, with a separation of 0.15" or 12 pc between 
the two isophotal peaks Lauer et al. (1996) and a black 
hole mass of M bh ~ 1O 9 M (Magorrian et al. 1998). The 
transition radius in this galaxy would be r t ~ 17pc; so 
again we find that the eccentric stellar disk lies within r t . 
From cquation(6) and using 12 pc for the disruption ra- 
dius we estimate that the progenitor core density would 

;er.ca/ harris/mwgc.dat 
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Fig. 1. — In this figure we outline the range of progenitor core radii and central densities that would allow the creation of an eccentric disk 
near a massive black hole. Dotted lines are for a black hole mass of M bh = 1O 7 M . The dashed lines are for M bh = 3 X 1O 7 M and the dot 
dashed lines for = 10 s Mq. Horizontal lines are set from the constraint < rt where we require that the disruption radius be within 

the region where the black hole dominates the gravitational potential (the sphere of influence). The allowed region for the progenitor core 
is that above the horizontal lines. Vertical lines are set from the constraint ro < rt, requiring that the core radius of the cluster is smaller 
than the region in which the black hole dominates the gravitational potential. The allowed region is that to the left of the vertical lines. The 
diagonal set of lines are set from the constraint < where we require that the mass of the active part of the disk ~ Po r o be less 
than the black hole mass. Parameter space for progenitor clusters that could result in the formation of an eccentric is limited to the upper 
leftmost area on the plot. Galactic globular clusters (Harris 1996) (assuming M/Ly = 5, solar units) are shown as triangles. The Gl cluster 
in M31 (labeled as Gl) is placed using parameters by Meylan et al. (2001). M31, M32 and M33 are placed on this plot with values derived 
by Lauer et al. (1998). In the case of M32 the core radius is only an upper limit. The Galactic center (denoted GC) is placed on this plot 
based on values summarized by Alexander (1999). Limits on the core of two additional power law galaxies, NGC 3377 and NGC 3115, with 
unresolved dense cores, are shown as diamonds with bars extending outward to show that these points are limits. The values for these two 
galaxies are based on values derived by Faber et al. (1997). 
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be po ~ 6 x 10 5 M Q pc~ 3 . The width of the peaks com- 
prising the disk (and hence the progenitor core radius) is 
estimated to be 2-4 pc. However, NGC 4486A is much 
more distant than M31, and consequently the eccentric 
disk mass is probably much higher than that of M31, prob- 
ably exceeding 10 8 M Q and this is too massive for a globu- 
lar cluster. However as in the case of M31's estimated pro- 
genitor, the central densities, core radius and total mass 
are reasonable for a galaxy core. 

2.4. Comparison to the properties of galaxy cores 

In equation(2) we have assumed an isothermal profile for 
the stellar density profile of the background galaxy. How- 
ever, as shown by multiple studies (e.g., Faber et al. 1997; 
Lauer ct al. 1998) galaxy cores at small radii are seldom 
fit with a surface brightness profile proportional to R^ 1 
(where R is the observed angular radius from the nucleus 
on the sky) corresponding to an isothermal density profile 
oc r~ 2 . Images from the HST have been part of a major 
effort to classify the nuclear stellar profiles in early-type 
galaxies, resulting in the classification of light profiles into 
two categories, galaxies with shallow inner cusps, denoted 
'core-type' profiles and galaxies with 'power-law' light pro- 
files (Lauer et al. 1995; Faber et al. 1997). Power-law 
type galaxy cores tend to have the steepest nuclear sur- 
face brightness profiles fJ,(R) oc i?~ 7 with 7 nearly equal 
to 1 in only the the most extreme cases. Most of the galax- 
ies studied by Faber et al. (1997) had shallower profiles. 
This implies that equation(3) somewhat underestimates 
the transition radius in most galaxies, particularly in the 
galaxies classified as core-type in which 7 is close to zero. 
We can regard our estimated radius (equation 3) as a lower 
limit on the transition radius for all but the galaxies with 
the steepest nuclear profiles. 

We now consider the problem of two merging galax- 
ies, both with more complex stellar density profiles and 
both with massive black holes. We can approximately de- 
scribe the primary as having a transition radius given by 
equation(3) (which is actually a lower limit as mentioned 
above). If the secondary has a power-law form for its den- 
sity profile then it will not completely disrupt at a par- 
ticular radius, like the non-singular isothermal sphere or 
King models. 

One component of the Nuker surface brightness pro- 
file (Faber et al. 1997) is the break radius, r^, where the 
slope of the surface brightness profile changes. For pro- 
files denoted core- type, within r^, the surface brightness 
(and so density) rises much less steeply with decreasing 
radius than outside it. Because the surface brightness for 
core- type galaxies only increases slightly within rb, we can 
associate the break radius and density estimated at this 
radius, with the core radius r and density po which we 
have used in the previous sections. This lets us deter- 
mine if the bulge or stellar component of the secondary 
galaxy will survive tidal truncation to within the transi- 
tion radius of the primary. A secondary with break radius 
exceeding its own transition radius rj, > r t (likely to be 
true for luminous galaxies) will almost completely disrupt 
when the tidal truncation radius reaches its break radius. 
The black hole of the secondary will then spiral in to the 
nucleus of the primary, leaving the stellar component of 
the secondary behind. 



However for bulges and fainter ellipticals, we expect 
rb < r t . For example, M32's profile is steeply rising well 
within its transition radius r t <~ 0.5pc and increases in 
density to its last measured point (Lauer et al. 1998). In 
this case the core would be continuously disrupted, though 
we expect a change in the surface brightness profile of the 
resulting disk set by the location at which the tidal trunca- 
tion radius of the secondary was equal to its own transition 
radius. The continuous stripping of the secondary in the 
presence of the secondary's black hole may prevent the 
formation of a lopsided disk. 

2.5. Tidal disruption of a power-law galaxy 

The radius of tidal disruption can be estimated by com- 
paring the mean density of the object to that of the ob- 
ject disrupting it (e.g., Sridhar & Tremaine 1992). At a 
distance r from the nucleus of the primary galaxy, the 
mean density resulting from the black hole is p(r) = 

Mbh,i/{^j-) where Mbh,i is the mass of the primary's 
black hole. We have assumed that r < r t .\ where r t .\ is 
the transition radius of the primary, so the contribution 
from the primary's bulge is negligible. 

For a secondary galaxy with a power-law law density 
profile (see Binney & Tremaine 1987 section 2.1e) 



(13) 



the mass integrated out to radius s is 
so the mean density within this radius 

*•>- <&*(!)" ■ < 15 > 

for radii s > rt,2 outside the transition radius of the sec- 
ondary. We set the mean density of the primary within 
r (due to the primary's massive black hole) equal to the 
that of secondary (within stripping radius s), p(s) — p(r), 
and find that 

3 M bh (3-a) 



(16) 



When the secondary is at a distance r from the nucleus of 
the primary, it will be stripped out to a distance s from 
its own nucleus, where s is related to r by the previous 
equation. 

The stripping radius, s, is the same size as the distance 
from the primary's nucleus, r, at a radius 

A 5 a ( 17 ) 

where Mbh,i is the mass of the primary's black hole. When 
s estimated from equation(16) exceeds r then the sec- 
ondary engulfs the nucleus of the primary. 

For a disrupted secondary galaxy core to produce an 
eccentric disk, we expect that the radius at which the sec- 
ondary engulfs the nucleus of the primary, r e , should be 
smaller than the transition radius of the primary (rt,i). 
In this case, the core of the secondary will survive intact 
within the transition of the secondary. A reasonable con- 
dition for formation of an eccentric disk via the disruption 
of a galaxy with a power-law density profile should be 

r e < r M . (18) 
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Using equations (17, 3) this is equivalent to 

Pa > ^a" a « a " 3 ^ (19) 

where k — 5.7 x 10 _4 pcM Q 1 ^ 2 . For an isothermal profile, 
the limit on the density at a = lpc 

/ M \ 1/2 
p a=lpc >5xlO«M 0P c-3 \^±L 



for a = 2. (20) 



2.6. Existing simulations of galaxy mergers 

With the exception of Merritt & Cruz (2002); Holley- 
Bockclmann & Richstone (2000), few simulations have 
been carried out with realistic stellar density profiles and 
central black holes. Merritt & Cruz (2002) presented INT- 
body simulations of two merging galaxies, the primary 
with a shallow central density profile and the secondary 
with a steeper profile at the nucleus. They did simulations 
with and without massive nuclear black holes to illustrate 
the difference in remnant profiles following the merger. 
The density profiles adopted for their simulations were 
generated from the family introduced by Dehnen (1993). 

To explore the role of the black holes on the disruption 
of the secondary in their simulations we estimate the tran- 
sition radii for each galaxy. At small radii, r < a, both 
galaxies were power-law in form 

where M was the mass of the galaxy and a was its scale 
length. Merritt & Cruz (2002) used ai = 1 for the pri- 
mary, and a 2 — 2 for the secondary (subscripts here refer 
to the primary and secondary galaxy respectively). We 
estimate the transition radius for each galaxy (sphere of 
influence) by calculating the radius at which the integrated 
stellar mass within equals that of the central black hole. 
For the primary galaxy 

f^iV /2 = 0.035 Ol 



V Mi 



(22) 



where we have integrated equation(21) and used the ratio 
M bh /M = 0.0012 (that adopted by Merritt & Cruz 2002). 
For the secondary galaxy we estimate 

' M bh i \ , N 

' ' (23) 



r t ,2 ~ a 2 



Mi 



= 0.0012a 2 . 



Because the transition radius of the secondary was so 
small, its black hole should not have affected the simula- 
tions (they were probably not well resolved on that scale). 

Both simulations produced similar remnants at radii 
larger than 0.1 times the primary's scale length. The 
change in the surface brightness profiles at this radius was 
probably due to the break radius (scale length at which 
the density profile changes slope) of the secondary. 

For the simulations containing both massive black holes, 
a change in radial density profile is seen in the remnant 
at approximately the primary's transition radius (see their 
figure 2d.) Using their profiles, we estimate that r e ~ 0.01 
times the scale length of the secondary, which is consis- 
tent with the survival of a fraction of the secondary's core 
to well within the transition radius of the primary. The 
constraints developed in the previous section were not vi- 
olated (r e < r ti i), so this type of simulation might have 



been capable of forming an eccentric disk in the remnant. 
If the simulation did not form one, then simultaneous dis- 
ruption of the densest part of the cluster might be required 
to form an eccentric disk. This would occur when the core 
has a well defined core radius (and nearly constant den- 
sity within), and may not occur when the secondary has a 
steeply increasing density profile all the way to its nucleus. 

In the simulations lacking the massive black holes (see 
their figure 2h), the secondary did not disrupt within its 
break or scale radius. This follows because the mean den- 
sity of the secondary significantly exceeds the mean density 
of the primary within their scale radii, consequently the 
primary could not have disrupted the secondary within the 
secondary's scale radius. 

The merger simulations of Hollcy-Bockelmann & Rich- 
stone (2000) were done with a quite massive primary black 
hole Mbh > 10 9 M Q and lacked a secondary black hole. 
From their initial density profiles we estimate that the size 
of the secondary at its disruption radius in their 10:1 sim- 
ulation was only about a third the disruption radius itself. 
This implies that a disk could be formed from the sec- 
ondary, following disruption. Hollcy-Bockelmann & Rich- 
stone (2000) reported that their simulations produced a 
spinning remnant disk, however it was was thick, consis- 
tent with our crude estimate of the secondary's largish 
size during disruption. Much of the secondary survived 
to within the transition radius of the pimary's black hole 
(about 30 pc). It is possible that a denser secondary (bet- 
ter satisfying the conditions in eqns 18-20) could result 
in a thinner stellar disk which might then be more likely 
to be lopsided. Alternatively, a secondary with a shallow 
core (that would would disrupt all at once, as simulated 
by Bekki 2000a) might be required for eccentric disk for- 
mation. 

3. DISCUSSION AND SUMMARY 

In this paper we have considered the possibility that a 
single disruption event could result in the formation of an 
eccentric stellar disk such as are found in M31 and possi- 
bly NGC 4486B. We explore the scenario of a galaxy core 
which is stripped as it spirals in toward a nuclear massive 
black hole until it reaches a critical radius, denoted here 
as rd, the disruption radius. In the simplest model of an 
isothermal sphere, the disruption radius is set by the clus- 
ter core or King radius and central density. We suggest 
that an eccentric disk cannot form unless the core of the 
galaxy or stellar cluster is disrupted within the sphere of 
influence of the massive black hole, referred to here as the 
transition radius, r t , and the core radius is smaller than 
the radius at which it disrupts. To succeed in making 
an eccentric disk, we find that the progenitor cluster core 
must be denser and more compact for less massive black 
holes. 

For the eccentric disk in M31, the progenitor cluster core 
density must be greater than 10 5 Mqpc -3 and core radius 
smaller than a pc. From the radius, width and mass of 
the most eccentric region of the disk we estimate that the 
density of the core was a few times lO 6 M0pc -3 and the 
core radius was of order a half pc. 

Massive extragalactic globular clusters such as Gl might 
be dense, compact and massive enough to be an eccentric 
disk progenitor, however they are probably too diffuse to 
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account for the large disk mass in M31 within a few parsecs 
of its black hole. 

Dense bulges such as found in M31 itself, and in lower 
luminosity ellipticals galaxies such as M32, are dense and 
compact enough that they could be similar to progenitors 
for M31's and NGC 4486B's eccentric disks. Nuclear star 
clusters, such as found in M33, may also be dense and 
compact enough to be progenitors. Though the mass of 
M33's cluster (a few million M©) is too small to be a pro- 
genitor for M31's eccentric disk, the nuclear stellar clus- 
ters observed by (Boker et al. 2002) in late-type galaxies 
have half light radii of order 5pc, range in their estimated 
masses from 10 6 — 10 8 M Q , and so could be progenitors 
for eccentric disks following disruption by a massive black 
hole. Because of the large mass of bulges and nuclear star 
clusters at small radii, it would be easier to account for the 
high masses of the two known eccentric disks by disrupt- 
ing them, than possible by disrupting a globular cluster. 
While the massive black hole in M32 might prevent a pro- 
genitor similar to M32 from forming an eccentric disk, nu- 
clear clusters such as found in M33 can lack massive black 
holes and so might provide better progenitor candidates. 
To date the high angular resolution of HST has resolved 
extremely high stellar densities or order 10 6 M Q pc -3 in 
only the nearest galaxy bulges. Until higher angular reso- 
lution observations are available, we will not know if such 
high density galaxy cores are common. 

The two galaxies with double nuclei, M31 and NGC 
4486B, exhibit only moderate color variations in their nu- 
clei (Lauer et al. 1996, 1998), a situation that could be a 
natural consequence of a scenario that involves the merg- 
ing of galaxy bulges (proposed here), and is more difficult 
to explain with a scenario that forms a younger disk in 
situ (an ingredient of the formation scenarios discussed 
by Bacon et al. 2001; Touma 2002; Taga 2002; Jacobs & 
Sellwood 1999). 



The scenario proposed here is based on a simple tidal 
disruption argument and can be most quickly tested with 
N-body simulations such as have been carried out by Mer- 
ritt & Cruz (2002); Bekki (2000a); Holley-Bockelmann & 
Richstonc (2000). The simulations of Bekki (2000a) es- 
tablished that the disruption of a cluster could result in 
the formation of an eccentric disk, and Merritt & Cruz 
(2002); Holley-Bockelmann & Richstone (2000) have car- 
ried out simulations based on realistic galaxy profiles and 
with massive black holes. Holley-Bockelmann & Richstone 
(2000) established that a spinning nuclear stellar disk can 
be a remnant. It remains to be seen whether the merger of 
two galaxy cores can result in the creation of an eccentric 
stellar disk. 

We suspect that the merger of a primary galaxy con- 
taining a very massive black hole with a secondary with 
a 'core-type' or shallow central surface brightness profile, 
and significantly lower mass black hole, would be most 
likely to form an eccentric disk with N-body simulations. 
For the core- type galaxies, the break radius and density 
at this radius provide an equivalent for the King or core 
radius and central density used in Figure 1 and so can be 
used to estimate the likelihood that the core can disrupt 
to form an eccentric disk. 

If the merger of galaxy bulges can result in the forma- 
tion of an eccentric disk, then their formation would be a 
natural consequence of hierarchical galaxy formation, and 
can also be used to probe the properties of the parents of 
galaxies which contain them. 

This work was initiated by discussions with Joel Green 
and Rob Gutermuth in the class Astronomy 552 at the 
University of Rochester during the fall of 2002. We thank 
Chien Peng, Ari Laor and Eric Emscllem for helpful dis- 
cussions and correspondence. 
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